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Abstract. The article presents a compact review of the analytical results (2002-2009) in 
the study of the system describing the motion of a top in two constant fields. The Liouville 
integrability of this system under certain condition of the Kowalevski type was established by 
A.G. Reyman and M.A. Semenov-Tian-Shansky. We present some geometrical foundations 
"^v^", of finding separations of variables. Two systems of local planar coordinates are introduced 

leading to separation of variables for two subsystems with two degrees of freedom in the 
dynamics of the generalized Kowalevski top. 
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1. Introduction 



Historically, the main goal in the study of a differential equation of mechanics was 
to integrate it, to find the explicit solution. The fact of the existence of principally non- 
integrable systems has shifted the accent to qualitative methods of investigation. However, 
during the last 20 years new algebraic approaches to construct or reveal the integrable 
systems were created. For Hamiltonian systems, many new cases of commutative and non- 
commutative integrability were found. A lot of them deal with the problem of the rotation 
of a rigid body about a fixed point and its generalizations to Lie groups that differ from 
5*0(3). Here lies one of the reasons for a new wave of mathematical interest in the integrable 
Hamiltonian systems. Note the books [U [2j [3l SJ |5] published lately and containing more than 
a thousand of references on the matter. In relatively recent studies, the integrable systems 
with three or more degrees of freedom were found not having continuous symmetries and 
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therefore not reducible as a whole to families of systems with two degrees of freedom (the 
so-called irreducible systems). One of the most ingenious and complicated example of the 
integrable irreducible system is the generalized Kowalevski gyrostat. It is a rigid body with 
a rotor satisfying the conditions of the Kowalevski type and placed in two independent 
constant fields, say, gravitational and magnetic. This case was found due to the efforts of 
several mathematicians. We note the works of O.I. Bogoyavlensky [6l[7J|8], H.M. Yehia [9|IT0]. 
L.N. Gavrilov I.V. Komarov [12]. The final and mostly general result was obtained by 
A.G.Reyman and M.A. Semenov-Tian-Shansky [131 [H]. Nevertheless, it appeared that no 
explicit solutions in the classical sense can be found for the general case. Up-to-date the 
only progress in studying the Reyman-Semenov-Tian-Shansky case is based on the idea of 
investigating critical subsystems. This row of publications starts from the work [T5], in 
which the Bogoyavlensky system generalizing the 1st Appelrot class is investigated without 
any connection to the general case, and the work [16] . in which the next generalization, 
of the 2nd and 3d Appelrot classes, was found. Later, the bifurcation diagrams of the 
partial cases and of the general case including the case of a gyrostat were built and classified 
in the space of parameters and with respect to the energy levels in the series of works 
[El [HI [191 EDI EH E21 E2]. The complete description was obtained for the critical set 
of the momentum mapping stratified according to various criteria [21]. For two critical 
subsystems generalizing the famous Appelrot classes separations of variables were found and 
the phase topology was investigated in [221 ESI [23 EH1 EH1 ED], all critical periodic motions 
were described analytically and integrated in elliptic functions for the cases of a top [211 EES] 
and of a gyrostat [321 EH]- From the integrability point of view, each of the critical subsystems 
found is as general as the whole Kowalevski case in the gravity field, i.e., generates a one- 
parametric family of Hamiltonian systems with two degrees of freedomj- These systems have 
all the features of the systems in which separation of variables is obtained algebraically. Here 
we call such systems algebraically solvable. We present though the geometrical ideas of how 
to construct such separations. 

A lot of publications deal with the notion of an algebraic complete integrable systems. 
This notion does not have a unique and exact definition. Usually it is associated with Hamil- 
tonian systems the topology of which is connected with the Jacobians and Prim varieties 
of algebraic surfaces. Algebraic complete integrable systems with known separation of vari- 
ables usually also have algebraic expressions of the initial phase variables in terms of the 
separated ones, although many contemporary papers, not like the classical works, do not 
pay much attention to this point. We think that a system with separated variables cannot 
be considered completely analytically integrated without explicit expressions for the phase 
variables. Separations of variables on algebraic curves generated by Lax representations, 
the existence of which is theoretically proved, do not provide any method to establish the 
connection with the initial variables. Here we use another approach — starting from the 
geometry of initial variables and trying to find such projections of the phase space which 
bring the integral manifolds to the most simple image. If it is done (and there are no general 
ways to do it), then the connection with the initial variables is kept from the beginning. In 
this paper we propose the projections and the coordinate systems that made it possible to 
find new algebraically and topologically interesting separations of variables. The examples 
given below show the valuable applications. 



In [33], all critical subsystems for the Kowalevski gyrostat in two constant fields are found having one 
extra parameter. 
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2. Geometrical approach to separation of variables 

2.1. Geometry of separated variables. The term we use here — algebraically solvable 
systems — is not a standard one. Let us give some explanations of the problem. Consider a 
system of differential equations on a submanifold V in a real vector space 

^ = X(x), (2.1) 

defining the evolution of the set of phase variables x and admitting separation of variables 

dsi 



(h v Pi(si,f), (2.2) 

where s« are so-called auxiliary or separated variables. We denote the vector of these variables 
by s. The functions Pi(s; f) are supposed to be polynomials in one variable s with coefficients 
depending on a set of arbitrary constants f , r is the "reduced time" connected with the real 
time t by the equation of the type dr/dt = p(s) > 0. Usually such separation has a 
mechanical origin, therefore equations (12. 2 p have the form of an energy integral. Moreover, 
we suppose that all phase variables Xj are expressed via the separated variables by rational 
functions with respect to the set of radicals 



Ria — \J Si — e a (e a G C) 

with coefficients smoothly depending on s. The set of numbers {e a } depending, of course, 
on the constants f includes also the roots of Pj. The classical examples show that the roots 
of Pi usually provide the whole set of {e a }. 

The vector f for the system (12. ip usually appears as a collection of arbitrary constants 
of some set of first integrals {J-'j}- The dependencies 

x = x(s;f), seAcc(f) (2.3) 

then are considered as the multi-valued parametric equations of an integral manifold JF{ = 
.F _1 (f) C V. Here Acc(f) is the region in the separated variables space filled with trajectories 
of (12. 2p for given f . In Russian, it is traditionally called by the term corresponding to "the 
region of the possibility of motions" [3"1"| [33] , which sounds not very elegant. A shorter 
term "accessible region" was introduced in English translations and completely reflects the 
situation. Notation of the region Acc(f) for the argument vector s in ( 12. 3 p is based on this 
term. 

Let the system ( 12. ip be a Hamiltonian complete integrable system with two degrees of 
freedom and 

J= : V ->• R 2 (2.4) 

be its integral mapping of compact character (the inverse image of a compact set is compact). 
Two degrees of freedom are chosen to simplify notation and because of the main examples 
belowS 

Definition 1. The set Im J 7 is called the admissible region {in the space of integral constants). 
A point f G R 2 is called an admissible point if and only if J-" _1 (f) ^ 0. 



Generally speaking, separation of variables is some mapping 

,2 



n:V^ RL S2) (2.5) 



In fact, here we do not use even the Hamiltonian type of the systems considered. 
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of the phase space to the plane of the auxiliary variables (si, s^) that takes the system (12. ip 
to the system (12.21) with i = 1,2. For any admissible f we obtain the mapping 



Definition 2. For a given f the set Acc(f) = 7Tf(J-f) is called an accessible region in the 
plane of the separated variables. 

It is clear that Acc(f) is a subset of 



and any point is included in Acc(f) with its total connected component in the set ( 12. 7p . 
Then Acc(f) is a collection of rectangular regions in the (si, S2)-plane; "infinite rectangles" 
(half-strips, strips or even quadrants) are also possible. The fact that the accessible region 
is rectangular for all integral constants, obviously, is the necessary condition of separation 
of variables. Possibly, supposing that the projection of the type (12. 6p is non-trivial in some 
sense or that the boundaries of accessible regions essentially depend on f in terms of the 
corresponding Jacobians, one can prove that this fact is also sufficient, but here we do 
not consider such a goal. We want to understand in which way to look for separations of 
variables. 

Note that separation of variables not always has the form of a global mapping (12. 5p . 
As far as the dynamics of a rigid body is concerned, separations not depending on integral 
constants are known for the Euler case (one independent variable can also be treated as 
a separation), the Goryachev-Chaplygin case and its generalization to a gyrostat given by 
L.N. Sretensky [36J, the Clebsch case with zero momentum constant (isomorphic to the 
Neumann problem). Already in the Kowalevski case [37], the change of variables leading 
to separation essentially uses the integral equations. However, having the mapping (12. 6p 
for all f one can formally substitute f = J-"(x) to obtain some mapping of the type (12. 5p . 
It is important to mention that practically in all cases separation of variables is not built 
straightforwardly in the form (I2.5p . As a rule, some intermediate mapping V — > R| is 
constructed first. Let us call this auxiliary plane of the variables ^1,^2 the picture plane. 
With the projection onto the picture plane, the separated variables are introduced as the 
functions of £i,£ 2 ,f (usually defined by some implicit equations). 

What is the geometry behind this construction? The mapping of the phase space to 
the picture plane takes an integral manifold to a region JFf bounded by some curves. If 
these curves could be included in a two-parametric family of curves on the plane R| and 
this family could serve as a local coordinate net, then in these coordinates the image of each 
integral manifold would be a rectangle. Therefore, the new coordinates satisfy the necessary 
condition of a separation and we obtain hope to build such a separation. At this moment 
it absolutely does not matter whether this coordinate net was built by a global change of 
variables or, as in the Kowalevski case, this net (the Zhukovsky net) depends in a very 
complicated way on the integral constants. 

Studying the projections of V onto the picture plane (or in more general case onto some 
picture manifold with dimension equal to the number of degrees of freedom) we must have 
in mind that V itself may not have a global coordinate system. For example, in the rigid 
body dynamics V = R 3 x 5*0(3) with the matrix group 5*0 (3) naturally embedded in R 9 by 
means of the so-called direction cosines. Any system of three coordinates either not covers 
SO (3) or has inevitable singularities (like the Euler angles or their modifications). It is one 



'(si,s 2 )- 



(2.6) 




(2.7) 
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of the reasons to consider V as a subset in a vector space V of sufficiently high dimension^! 
presenting it as the level surface of some mapping 

Y . ~y y j^dim V— dimP 

In the rigid body dynamics this mapping is formed by the so-called geometric integrals, 
and in the theory of the Euler equations on Lie algebras this role is given to the Casimir 
functions. 

Combining (12 .4p with T we obtain each integral manifold as a level surface of the new, 
enhanced, integral mapping Q : V — > Z, where V and Z are real vector spaces and J-f 
coincides with (? z =(? _1 (z), z G Z. By an appropriate choice of coordinates we can represent 
V = XxY, where X is the picture plane, and replace the projection of V onto the picture 
plane with the projection of V onto the first multiplier px '■ XxY X . This construction 
is easily generalized to the case when X, Y, Z are smooth manifolds. In axially symmetric 
problems of the rigid body dynamics the manifold X is often identified with the so-called 
Poisson sphere S 2 . Then instead of Acc(f) we have a similar region 

Acc(z) = lmp x \g z = p x {Q" 1 {2,)). 

To construct coordinate nets on the picture plane we need to investigate the boundaries of 
accessible regions, not only in the topological sense, but the curves that in geometry and in 
the wave fronts theory are called visible contours. 

Definition 3. [35]. The generalized boundary of an accessible region Q x is the set of critical 
values of the mapping 

Px\g x :Gz^X. 

To find generalized boundaries one needs, in particular, to solve the equation Q = z. 
The following result [35] helps to avoid this sometimes impossible procedure. 

Proposition 1. Let v 6 V. Denote by L(v) the restriction of the operator 

T V G : XxY Z 

onto the subspace Y . The image of the point v e Q z belongs to the generalized boundary of 
the accessible region Acc(z) if and only if rankL(v) < dimZ. 

Naturally, this statement is of practical use only if dimQ z ^ dimX, which is equivalent 
to dimY" ^ dimZ. In global coordinates on V and Z, the condition for rankL(v) is easy to 
write in terms of the corresponding Jacobi matrix. 

2.2. Integrable problems of the rigid body dynamics. We tend to apply the above 
ideas to the system of equations describing the rotation of a rigid body about a fixed point 
in the linear potential force field 

T du> _ 
I—— = la; x oj + r x x ck + r 2 x p, 

dS d(3 a ( 2 - 8 ) 

— = axw, — = pxu. 
dt dt 

Here oj = (u>i, <jj 2 , uj 3 ), a = (a±, a 2 , « 3 ), (3 = (j3x, 02, ^3) are the phase variables, the constant 
vectors r 1; r 2 G R 3 and the diagonal constant matrix I = diag{A, B, C} are the physical 
parameters. The geometric integrals forming the corresponding mapping T : R 9 — > R 3 are 
2 ,f3 2 ,a-f3. These are the Casimir functions for the Poisson brackets on R 9 bringing the 



The main reason is, of course, to have a possibility to introduce algebraic and rational functions. 
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equations (12. 8p to the Hamiltonian form [8]. If at, f3 (the intensity vectors of the forces) are 
linearly independent, the common level of the geometric integrals is a smooth 6-dimensional 
manifold V endowed with a symplectic structure. 

Put I = diag{2,2, 1}, n = (1,0,0), r 2 = (0,1,0). Then equations (J22D describe the 
Kowalevski type top in a double force field. The first partial case of integrability (with 
restrictions on the phase variables) was found in (6J [Tl [8] . In the special case of Yehia [9] 
(with restrictions on the force fields \ot\ = \f3\,a-f3 = 0) the system has a symmetry group 
and becomes reducible to a family of integrable systems with two degrees of freedom. The 
complete integrability of this problem in general was established in [T3| [TJ]. More detailed 
presentation for the case of the Kowalevski gyrostat is given in |38j . 

Definition 4. Let ct, (3 G R 3 . The following constants 

p= sJa 2 + (3 2 > 0, r= \j{ot 2 - (3 2 ) 2 + A(ot-(3) 2 ^0 (2.9) 

are called the invariants of the pair (a.,f3). 

Obviously, p ^ r and p = r if and only if ckx/3 = 0. It is the classical Kowalevski case. 
If r = we obtain the reducible Yehia case. 
For 6 G R denote 

e= (cos9 -sin0\ 0=(4|fV (2.10) 
y smO cos J V 1 J y ' 

Proposition 2. [T5] The linear automorphism o/R 9 

(£) =e ()) §T ' u '= §w - < 2 - n > 

preserves the equations (12.81) . An orbit of the restriction of this automorphism onto the space 
R 6 = {(ct,f3)} consists of all pairs having the same invariants. Each such orbit contains an 
orthogonal pair ot f3 = 0. 

It follows that without loss of generality for irreducible problems we can consider the 
phase space V C R 9 as given by the equations 

a 2 = a 2 , (3 2 = b 2 , a-(3 = (a > b > 0). (2.12) 

The proposition holds for an axially symmetric gyrostat in the same force field [23J. 

The integration of the whole generalized Kowalevski system has not been fulfilled. Sep- 
aration of variables of the algebraic type is obtained only in the following three cases: 

• the subsystem with one degree of freedom consisting of the families of critical periodic 
motions in the first critical subsystem with two degrees of freedom; the subsystem 
itself was found by O.I. Bogoyavlensky in the works [SI E], the existence of periodic 
motions was proved in [15], integration presented in [3"T] ; 

• the second critical subsystem with two degrees of freedom found in [T6j , separation 
of variables is pointed out in [251 EE] ; 

• the third critical subsystem found in [18], separation of variables is given in [28 | [29] , 
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For the generalized Kowalevski top, let us present the system ( 12. 8 j) in the scalar form 

(2.13) 



UJ 2 UJ 3 + @3, 


2lo 2 = 


-uj 1 uj 3 - 


~ «3, 


tb 3 = a 2 


Oil = «2^3 


- "3^2, 


$1 = 


ft^ 3 


- ftw 2 , 


a 2 = a 3 u 1 


- "1^3, 


k = 


ft^l 


- P1W3, 


a 3 = a x u 2 


- a 2 u ± , 


ft = 


01^2 


- ft^i- 



(2.15) 



Its first integrals £1X6 £IS follows P [Lj 

if = u 2 + u 2 + ^uj- K + ft), 

if = (w 2 -ui% + ax- ft) 2 + (2wiw 2 + a 2 + ft) 2 , 

1 1 ( 2 - 14 ) 

G = {aiUi + a 2 ui 2 + -a 3 u 3 ) 2 + (ftwi + ftw 2 + -ftw 3 ) 2 

1 

+ ^3(71^1 + 72W2 + -73^ 3 ) - «i6 2 - fta 2 . 

Here 7, stand for the components of the vector ax/3. We denote the corresponding constants 
by h,k,g. 

To make the description of critical subsystems compact, we use the complex change of 
variables [18] generalizing that of S.V. Kowalevskfl, 

aci = (an - ft) + i(«2 + A), £2 = ("i - ft) ~ i(«2 + ft), 
?/i = (aii + ft) + i(a 2 - ft), V2 = ("i + ft) - i(a 2 - ft), 
21 = a 3 + ift, 22 = "3 - ift, 
w 1 = u 1 + ioo 2 , w 2 =u 1 - ioo 2 , w 3 = u 3 , 

where i 2 = — 1. It brings the first integrals to the form 

H = \ W l + W 1 W 2 ~ + 2/2), 

K = (w\ + xi)(wl + x 2 ), 

G =^(p 2 - xxx^wl + -(x 2 z 1 Wi + x l z 2 w 2 )w 3 (2.16) 

+ ^{x 2 wi + Vxw 2 )(y 2 w l + xiw 2 ) - ^P 2 {yi + 2/2) 

+ ^r 2 (xi +x 2 ). 

Here according to $ZM), (l2~12|) . 

pi = a 2 + b 2 , r 2 = a 2 -b 2 . (2.17) 
Equations (12.121) take the form 

z 2 + xiy 2 = r 2 , z 2 + x 2 y 1 = r 2 , (2.18) 
x 1 x 2 + y 1 y 2 + 2z 1 z 2 = 2p 2 . (2.19) 

The first critical subsystem A4i C V is defined by the equations 

w\ + Xi = 0, u>2 + x 2 = 0. 

On Mi we have 

k = (2.20) 



4 These combinations of the real variables are in fact found in the Lax representation [T3] . 



8 MIKHAIL P. KHARLAMOV, ALEXANDER Y. SAVUSHKIN 

and we take H and the Bogoyavlensky integral 

F = wiw 2 w 3 + z 2 wi + z 1 w 2 

as a pair of almost everywhere independent integrals. 

The second critical subsystem A^2 C V is defined by the equations 

xix 2 w 3 - (x 2 z 1 w 1 + x 1 z 2 w 2 ) = 0, —(wl + xi)-—(w 2 + x 2 ) = 0. (2.21) 

x 1 x 2 

On Ai 2 , we choose the following pair of almost everywhere independent integrals: 

M = ^- 2 [^(wl + x 1 ) + ^(w 2 2 + x 2 )], 

I Xl 002 (2.22) 

L = (wiw 2 + x x x 2 + z x z 2 M). 

y/XiX 2 

The constants of the general integrals (I2.14p satisfy the equation 

(p 2 h - 2g 2 ) 2 - r 4 k = 0. (2.23) 
The third critical subsystem Ai 3 C V is described as 

W 2 X\ + lOi2/ 2 + ^3^1 WiX 2 + 10 2 2/1 + ^3^2 



0. 



Wi w 2 

v2 „„ Jl 



(W 2 ^l + Wl2 2 )w3 + 



+Xiiul + x 2 w\ 



w 2 zi w x z 2 , . 

L iOi u> 2 (2.24) 

wix\Z\ w\x 2 z 2 
w 3 + - 



w 1 w 2 
+x 1 z 2 w 2 + x^iwx + (wiz 2 - w 2 z 1 ){yi - y 2 ) = 0. 

Here we take as a pair of almost everywhere independent integrals 

s _ 1 , Wjh + x l w 2 + ^1^3 ^2^1 + V\W 2 + 2: 2 W 3 x 

4 ^ w\ w 2 ' ' 

^ ! r / s \ i (2-25) 

T = ^FH^wi + yiiw 2 + z 2 w 3 ) + w 2 {y 2 wi + xiu> 2 + ^1^3)] 

+ XiX 2 + z x z 2 . 

The constants of (I2.14p on each manifold {S = s, T = r} satisfy the following equations 
p 2 — r , r 2 — 2» 2 r + r 4 p 4 — r 4 1 . 9 

ft = V + s - *= — h — + T - ° = — + 2 ip - T)a - (2 - 26) 

As b tends to zero the critical subsystems turn into the famous classes of "especially 
remarkable" motions found by G.G. Appelrot |39j in the Kowalevski case. From the point of 
view of the phase topology of the integrable Hamiltonian system these critical subsystems 
form the critical set of the momentum mapping and equations (12.201) . (I2.23P and (I2.26P define 
the bifurcation diagram 



2.3. Two systems of local coordinates. 

2.3.1. The coordinates si,s 2 . We follow the advice of C.G.J. Jacobi [40J to find first a con- 
venient change of variables and then search for the problems in which it may be successfully 
applied. 

Consider a number r > 0. On the plane (x, z) introduce two one-parametric families of 
circles 

x 2 + z 2 + r 2 = 2six, x 2 + z 2 - r 2 = 2s 2 x. (2.27) 
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Except for the points of the axes Ox, Oz these families provide the net of curvilinear coor- 
dinates (s%, s 2 ) on the (x, z)-plane. In Fig. [T]we show: (a) the curves Si = const for Si ^ r; 
(b) the curves s 2 = const for s 2 in some interval about 0; (c) joint net (si, s 2 ) scaled 1:2 for 
better view. For the differentials we have 

/£j2 ^2 ^2 ^ ^2, ^,2 j ^2 ^ 

dsi = dx H — dz, ds 2 = 7. dx H — dz. (2.28) 

2x 2 x 2x 2 x 




(c) 

Figure 1. The coordinate net (si,s 2 ). 

We now establish the connection with the rigid body configuration space (I2.12p and find 

the image of it on the (si, S2)-plane. To this end, choose the parameters p,r according to 
f)2.17p . Introduce the variables follows 

x 2 = xix 2 , y 2 = yiy 2 , z 2 = z x z 2 . (2.29) 

Then due to (12.191) we obtain the ellipsoid 

x 2 + y 2 + 2z 2 = 2p 2 . 

From (12.191) express 

y lV2 = 2p 2 -x 2 - 2z 2 (2.30) 

and represent (12.181) in the form 

(zi + z 2 f = 2r 2 - ( Xl y 2 + x 2 y t ) + 2z 2 , ^ ^ 

{z 1 - z 2 f = 2r 2 - ( Xl y 2 + x 2 y x ) - 2z 2 . 1 ' ' 

Elimination of the variables z±,z 2 and the product y\y 2 in (I2.18p . (I2.19p . (I2.30p gives 

r 2 [x x y 2 + x 2 yi) = r 4 + 2p 2 x 2 - {x 2 + z 2 ) 2 . (2.32) 
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Denote 

$±(x, z) = (x 2 + z 2 ± r 2 ) 2 - 2(p 2 ± r 2 )x 2 . 
Then (I2.3ip . ( I2.32p yield the following expressions: 

r 2 (z! + z 2 ) 2 = $+(x, z), r 2 (z x - z 2 ) 2 = $-(x, z) 

Therefore, the region for x, z is defined by the inequalities 

$ + (x,^)^0, $_(x,z)^0. 

Define s±, s 2 as in (12.271) with the chosen r. Then 



(2.33) 
(2.34) 



si 



x 2 + z 2 + r 2 



S2 



x 2 + z 2 — r 2 



2x 2x 

The conditions ( I2.34p take the form 

2 \ 2 2 ^ 7,2 

This rectangular area (in the above mentioned generalized sense) is the image of the phase 
space V under the projection (u>, a, (3) !->■ (a, /3) (si, s 2 ). The image of V on the (x, 2;)- 
plane with the corresponding part of the (s\, S2)-net is shown in Fig. |2]for the first quadrant. 



(2.35) 
(2.36) 




Figure 2. The image of the phase space. 



2.3.2. The coordinates ti, t 2 . Consider the plane with coordinates (£, x) and the second order 
curve 

£ 2 x 2 

^ + — = 1, (2.37) 
a t 

where r, a are not simultaneously negative. Let (£, x) be an arbitrary point of the region 

r£ 2 + ax 2 > to. (2.38) 

Then from this point we can draw exactly two tangent lines to the curve (I2.37p . In other 
words, the system of tangent lines to the curve (12.371) provides, in the region (I2.38p . the net 
of straight lines such that each point belongs to exactly two lines of this net (see Fig. [3]). 
This fact is well-known in geometry. 

To describe this net, introduce the variables ti,t 2 as the roots of the quadratic equation 

o 2t£ t£ 2 + ox 2 , 
t 2 ^t+— 5- = 0. 2.39 

T — X z T — X z 

Its discriminant is positive in the region f !2.38j> . hence the following expressions are real, 

tl = H±i£ t 2 = T ±Z^. (2.40) 
r — x z t — x z 
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Figure 3. The net of tangent lines to a second order curve. 



Here fi = a/ r£ 2 + ax 2 — to G R. 

Let us also introduce the connection between two coordinate systems £2) and (si, s-i). 
It will be useful below. Let p, r, a, b satisfy (12.171) . In (12.371) we put 

a = t 2 - 2p 2 r + r 4 

and let the coordinates £, x depend on (s\, S2) 

£=£i±iV-r, x = -^-. (2.41) 

Si — S2 Si — S2 

With respect to the variables (I2.29P we get 

£ = ^ + 2 2 - r. (2.42) 
Then the boundary lines of the region (12.361) 

si = ±a, s 2 = ±6 (2.43) 
appear to be coordinate lines in the net (ti,t 2 ): 

ti = -r±r 2 , t 2 = -r=Fr 2 . 

3. Example of a geometry based elliptic separation 

Here we present the result announced in [25, 26J. Complete investigation including the 
study of the phase topology can be found in [27] . 

Consider the second critical subsystem and fix the constants of the first integrals ( 12. 22ft : 

M = m, L = £. (3.1) 

Eliminate from (I2.2ip . ( 12 . 18j) . (I2.19P and (13. ip the variables w±, W2, W3, yi, y<i and the sum 
z\ + Z2- The result is one equation containing X\,X2 and the product z\Z2'- 

1 r -1 2 

m{x\X2 + Z1Z2) + v/ r A m 2 — r 2 ra{x\ + x 2 ) + XiX 2 = £ 2 . (3.2) 

x x x 2 L J 

It defines a two-dimensional surface in R 3 (3mxi, D\z X\,z). Introducing the polar angle by 
putting x\ + x 2 = 2xcos0 we obtain the equation on the cylinder R 2 (x, z)xS 1 (cf>mod27i): 

m[x 2 + z 2 ) + a/ r 4 m 2 — 2r 2 mx cos + x 2 ± £x = 0. 
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According to Proposition (TJ the generalized boundary of the projection onto the (x, z) -plane 
of a set given by an equation of the type f(x, z, 0) = without any restrictions on is 
defined by the condition df/d<f) = 0. In our case since the main equation contains only 
cos0, this condition immediately gives sin0 = 0. Then on the generalized boundary 

m[x 2 + z 2 ) + \r 2 m ± x\ ± tx = 0. 

Hence 

x 2 + z 2 ± r 2 ± ^-^x = 0. 

m 

Recalling the definition ( I2.27P we see that for any combination of signs these curves on 
the (x, z)-plane become the coordinate lines of the (si,S2)-net. In particular, all integral 
manifolds (13. ip have rectangular images on the (si, s 2 )-plane (possibly unbounded if the 
value x = equivalent to si = oo is admissible). Let us show that this fact indeed leads to 
separation of variables. 

Write the equations (I2.13P for the variables (12. 15ft : 

x[ = —X1W3 + ZiWx, x' 2 = x 2 w 3 - z 2 w 2 , 



y[ = -yiw 3 + z 2 w u y' 2 = y 2 w 3 - z x w 2 , 

2z[ = x\w 2 - y 2 wi, 2z' 2 = —x 2 w\ + yiw 2 , 

2w[ = -wi(w 3 - A) - z h 
2w' 2 = w 2 (w 3 - A) + z 2 , 
2w' 3 = y 2 -yx. 



(3.3) 



Here prime stands for d/d{it). Calculating the derivatives by virtue of these equations we 
get 

:{xiz 2 w 2 - x 2 z 1 w 1 ) 1 



x' 2(xix 2 ) 3 / 2 

(— )' = 2jx~x~y^^ XlZl ^ 2 + V2Xl ^ W2 ~ X2Z2 ( Z ± + yi x 2)w\]- 

Taking into account (13.11) . we obtain from (|2~3 



(3.4) 



So 



4^ l + ^k/^ 2 -\/^ l] * 



.2 



(3.5) 



Denote 

and express from (12.211) 



i?i = \/ r 2 m — Xi, R 2 = \/ r 2 m — x 2 



X\ x 2 

— w 2 = Ri, a — wi = R 2 . (3.6) 
x 2 V x i 



Equation ( 13. 2 p takes the form 



m[x 2 + z 2 ) + R Y R 2 = Ix. (3.7) 
It follows from the definition of Ri, R 2 that 

R 2 + R 2 = 2r 2 m - (x 1 + x 2 ), 



R\R 2 = r A m 2 — 2r 2 m(x 1 + x 2 ) + x 2 



(3. 
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Therefore from (13. 7p we obtain 



R1R2 = £x — m(x 2 + z 2 ), 
1 



R\ + Rl = ^{l £x - m ( x2 + z2 )f ~ x2 } + r<lm - 

Using the expressions of x, x 2 + z 2 via s\, s 2 , we find 

2 

(R x ± R 2 ) 2 = -{{I 2 - 1) - 2£m[( Sl + s 2 ) 

m{si - s 2 y 

- s 2 )] + 2m 2 [(s 2 + s 2) -p (s 2 _ s 2 )]}> (3 _ 9) 
Introduce the following notation: 



*(si, s 2 ) = 4msis 2 - 2£(si + s 2 ) + — (£ 2 - 1), 

m 

$( s ) = *( s , s) = 4ms 2 - 4£s + — {£ 2 - 1). 



Hence 



(zi + z 2 Y = ^(sf - a 2 ), ( Zl - z 2 ) 2 = -(s 2 - b 2 ). 



2r rz : 2r 



z 1 + z 2 = J si -a 2 , zi-z 2 = V-sl-fe 2 , (3.12) 

Si - S 2 V si - s 2 



(3.10) 



m 

Then from ( l3T9|) . 

Ri + R2 = V*(s 2 ), Ri-R 2 = —t—VHsi) (3.11) 

Si — s 2 Si — s 2 

(the radicals may have arbitrary signs). From (12.331) we get 

4r 2 , o ON , NO 4r 2 



where the signs of the radicals, similar to ( 13. lip , are arbitrary. 

Substitute (13. 6p . (I3.10p . (13.121) in (13.51) and return to the real time derivative to obtain 
the equations describing the time evolution of the auxiliary variables si, s 2 : 

f = iy^V^0, ^ = (313) 

In these equations the variables are separated and the equations themselves are easily inte- 
grated in Jacobi's elliptic functions. 

To accomplish the separation, it is necessary to present the expressions of all phase 
variables in terms of the separated ones. The values z±,z 2 are found from (I3.12p . From 

(E2ZD 

r 2 2 2 2 Si + s 2 

x = , x + z = r , 

Si — S 2 Si — s 2 

then xi,x 2 are found from (13.81) and (13.111) . The equations (I2.18P serve to define yi,y 2 , and 
w\ , w 2 are expressed from ( 13. 6 p . Finally, the variable W3 is defined from the first equation 
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(I2.2ip . Thus, for the complex phase variables we have 



and 



Denote 



Then 



~,2 



„2 



r 



x 2 



2( Sl -s 2 ) 2 



n (2s l32 - p 2 ) -2y/(sj- a 2 )(sl-t?) 
Hi — 2 / 

*(S1,S 2 ) - V $ ( S l) $ ( S 2) 



o (2 glg2 -p 2 ) + 2 v /(^-a 2 )( g i-^) 
1/2 = 2 . 

*(S1,S 2 ) + V $ ( S l) $ ( S 2) 



Si - s 2 v v 
Sl — s 2 v v 



Ml 



*(si,s 2 ) - v/$(si)$(s 2 )' 



u> 2 = r 



*(si,s 2 ) + v /$(s 1 )$(s 2 )' 

1 



si — s 2 



Si = \ s\-a 1 1 Vl = J - $( Sl ), S 2 = \b 2 -sl (p 2 = J$(s 2 ). 



a i = 77? T?[( s i s 2 - a + SiS^i^], 



2(si - s 2 ) 



«2 = - - ^ 2 [(sis 2 - a 2 )y?i^2 - S , i5 , 2 *], 
^ = ~ Ws~i - s 2 ) 2 ^ lS2 ~ ^ViV^ - S'i5 , 2*], 

& = 2(si ^ g2)2 [(sis 2 - 6 2 )^ + S 1 S 2 <p 1 <p 2 ], (3.14) 

«3 = Si, /3 3 = S2, 

Si - S 2 Si - s 2 

^1 = 2(3 - g ) ^ ~ 2mg i)y2, ^2 = ^ - - r (f-2ms 2 )yi, 

u 3 = — - — (S 2 <fi - Si<f 2 ). 
Sl - s 2 

Together with the differential equations (I3.13P this gives the complete analytical solution of 
the problem. 

4. Example of a hyperelliptic separation 

4.1. The geometry of accessible regions. Let us consider the dynamical system defined 
by the equations (I2.13P on the invariant almost everywhere four- dimensional manifold .M3 
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cut in R 9 by five equations (I2.18p . (I2.19p . (I2.24p . This system is Hamiltonian with two de- 
grees of freedom except for the zero-measure set of points where either the equations (12.241) 
fail to be independent, or the induced symplectic structure degenerates [28]. The possibility 
of separation of variables was predicted in [2T] basing on the above shown geometrical spec- 
ulations, but it took a long time to obtain the final analytical expressions, complex first [28] 
and finally real [29J. 

The functions S, T defined by (I2.25P provide the complete set of integrals in involution 
on In particular, the equations of the integral manifold 

{(eV:H(0 = h,K(0 = k,G(() = g} 

with the functions (I2.16P are replaced with the invariant relations (12.24(1 and the system of 
equations 

S = s , T = r. (4.1) 
Using the variables ( I2.15P we transform the system (12.24)1 . ( 14. ip to the form 

(y 2 + 2s)wi + xiw 2 + ziw 3 = 0, 

x 2 w x + (y 1 + 2s)w 2 + z 2 w 3 = 0, , ^ 

x%Z\W\ + x\z 2 w 2 + (r - x\_x 2 )wz = 0, 
2swiw 2 - (xix 2 + z x z 2 ) + r = 0. 

Using the variables (12.15)) with the corresponding dynamics (13.31) note that despite the 
complexity of these variables their space is 9-dimensional since the pairs (xi,x 2 ), (yi,y 2 ) 
and (zi,z 2 ) are complex conjugate and is in fact real. Seven relations (14. 2p . (12 . 18)1 and 
(12.191) define in this space the integral manifold Ts^. If the latter set does not have any points 
of dependence of the functions S, T, then it will consist of two-dimensional tori bearing the 
quasi-periodic motions. 

From what follows, we exclude the four points of equilibria on V. At non-trivial motions 
we have u) ^ for almost all time moments. Then the determinant of the first three equations 
(14. 2 p in Wi (i = 1, 2, 3) is identically zero. Eliminate in this condition the variables zf, z\ and 
the product y x y 2 with (I2TT8D . (12^5(1 and Q2J50]) to obtain 

2s[{r 2 x 1 - ryi) + (r 2 x 2 - ry 2 )\ = -r 2 {xYy 2 + x 2 yi) (A „x 

+2[2s 2 (r - x 2 ) + p 2 {r + x 2 ) - r[x 2 + z 2 )]. [ ^ 6) 

On the other hand, it follows from QZgJ) , (I2^D]) that 

{r 2 Xi — ry 1 )(r 2 x 2 — ry 2 ) = r A x 2 + r(2p 2 — x 2 — 2z 2 ) — r 2 r(xiy 2 + x 2 yi). (4.4) 

Denote 

a = t 2 - 2p 2 r + r 4 , x = ^ 0. (4.5) 
We have, from the second equation (12.26)1 . 

As 2 x 2 = cr + 4s 2 r. 
Introduce the complex conjugate pair 

/ii = r 2 xi - ryi, fi 2 =r 2 x 2 -ry 2 (4.6) 

and denote for brevity 

f = x 2 + z 2 - r. (4.7) 
Eliminate in flOJ), (1Q(1 the term x x y 2 + x 2 y x by (l2~32l : 

2s( y u 1 +/i 2 ) = £ 2 -As 2 (x 2 -t) -a, /i ly u 2 = r£ 2 + ax 2 - to. (4.8) 
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These two equations in four variables describe the integral manifold J- 8)T in the space 
(x, £, fjLi, /i 2 ) (a*2 — W)- According to Proposition [fl at the points covering the general- 
ized boundary of the projection of J~ S)T onto the (x, £)-plane the Jacobian of the system with 
respect to the variables //i,/i 2 equals zero, or in view of ( 14. 5p . 

(f - 2sx - 2s X )(£ - 2sx + 2s X )(£ + 2sx - 2s X )^ + 2sx + 2s X ) = 0. 
So, the straight lines 

£±2sx±2sx = (4.9) 

together with (12.43P bound the projection of the integral manifold onto the (x,£)-plane. 
Comparing ( 14. 7p with (12.42p recall that the boundary lines of the region (I2.36P after the 
change (I2.35P become tangent to the second order curve 

r£ 2 + ax 2 - ra = 0. (4.10) 

The resultant in £ of any of the equation (14. 9 p with ( I4.10p depending on the chosen signs 
evaluates to 4s 2 (r ± x x ) 2 ■ Hence all straight lines (14. 9p are also tangent to the curve ( I4.10p . 
namely, at the points with x = =Ft/x. 

Choose Ai = y/2sfii + a, A2 = y/2sfi 2 + cr to be complex conjugate. Then the system 
( 14. 8 p simplifies 

(Ai + A 2 ) 2 = E+(x, z), (Ai - A 2 ) 2 = H_(x, z). 

Here 

S±(z, z) = (x 2 + z 2 - t ± 2sx) 2 - 4sV, 

and the solution exists if and only if 

E + {x,z)^0, E_(x,z)^0. (4.11) 

The inequalities (I2.34p . (14. lip describe the corresponding accessible region Acc(s,r) on the 
(x, 2;)-plane. We now present this region in the plane of the variables (si, S2). Express 



T 



X + Z - T = [Si + S 2 ^(Si - S 2 )]x, 

E+(x, z) = x 2 A + A_, EJx, z) = x 2 M,M_ 



(4.12) 



where 



A±(si, s 2 ) = si + s 2 - r f SX ( s i ~ s 2) ± 2s, 

t + 2sx 

M±(si, s 2 ) = si + s 2 -g — (si - s 2 ) ± 2s. 

From ( fOTjl . (I47T2]) we have 

A + ( Sl , s 2 )A_(s ls s 2 ) ^ 0, M+( Slj s 2 )M_(si, s 2 ) ^ 0. (4.13) 

The lines Aj. = 0, M± = form a parallelogram. The solutions of (14.131) fill two 
half-strips attached to its sides belonging to the lines A± = 0. Moreover, the solutions of 
(12.361) fill two horizontal half-strips attached to the sides s 2 = ±6 of the rectangle with 
vertices S\ = ±a, s 2 = ±b. The intersection of these sets give the accessible region for the 
variables Si, s 2 (e.g. see Fig. S]for a — 1, b — 0.4, r = 1.2, s = —0.6). With respect to these 
variables all boundary lines are also tangent to the second order curve, which is the image 
of the curve (I4.10p under the change (12.411) . 

The presented geometry on the planes ( £) and (si, s 2 ) shows that it is possible 

to separate the variables by use of the second coordinate system (12.401) . where the notation 
have the sense given in this section. 
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Figure 4. Example of an accessible region. 



4.2. Parametric equations of integral manifolds. At the two-dimensional surface ( 14.1 p 
the phase variables (12.15P can be expressed in terms of two. At this moment it is convenient 
to choose for that the variables x, £. Then from the equations (14. 8 p we have 

_ 2st) + (/i 2 - 2sr) = |- - 2s(x 2 + x 2 ), 
(//! - 2sr)(// 2 - 2sr) = As 2 X 2 x 2 . 



Hence 



/ii = 2sr + ^( v / ^+ v 7 ^) 2 , /i 2 = 2sr + — 
8s 8s 



(4.14) 



'1 — V ^2; 

8S SS 

Here we denote 

= e - 4s 2 (x + X )\ M*, = ^- 4s 2 (x - xf- 

From (12.301) we have 

ym = 2 P 2 -2(e + r) + x 2 . 

Consider this equation together with the definition X1X2 = x 2 and the equations (14. 6p . In 
view of (14.141) we get 



Xi 



x 2 



2s4r 4 (x 2 - t) + r(V$i 



r 2 16s 2 r + ( v /^ 1 - 
2s 4r 4 (x 2 - r) + r( y 7 ^ - y^) 2 
r 2 16s 2 r + lyMy + a/^) 2 



(4.15) 



2/1 = 2s 
1/2 = 2s 



4[2r£ - t(x 2 - r) + a] - (ygT - y 7 ^) 2 

16s 2 r+ (y 7 ^- y^) 2 
4[2r£ - r(x 2 - r) + a] - (y^ + y 7 ^) 2 



(4.16) 



16s 2 r + (V*i + y 7 ^ 



Here the polynomials 



$!(a:,f) = (£ + r + r 2 ) 2 -2(p 2 + r 2 )x 2 , 
$ 2 (x,0 = (£ + r - r 2 ) 2 - 2(p 2 - r 2 )x 2 



(4.17) 
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are obtained from $± after substituting z 2 with its value from ( 14. 7p . In the sequel such 
substitution is supposed by default. From (12.3 ip we then get 

zi = — {y/^ + yf¥ 2 ), z 2 = — (^1^ ~ V^2). (4.18) 

As a result we have the expressions for all phase variables from the so-called configuration 
group, i.e., determining the orientation of the body. The radicals in ( I4.15p . (I4.16p . (I4.18P 
have arbitrary signs. 

To find the variables W{, which are in charge of the angular velocity components, we use 
the expressions (12.161) for the general integrals K , H and the relations (I2.26P , which define 
the constants k, h in terms of s, r. In view of (14. 5 p we have 

K= {w 2 l +x 1 )(w 2 2 + x 2 ) = X \ (4.19) 
H = -wl + w l w 2 --{y l + y 2 ) = s + ^-^. (4.20) 



From (I4.19p . (I2.25P we obtain the system 

x 2 w\ + x x w\ = --^[£ 2 + 4s 2 (y 2 + x 2 )}, wxw 2 = J-. 

As z 2s 

Hence 



W ^4^ (V ^ + ^ "* = 5^^-^ (4 ' 2i; 



where 

= (£ - 2sx) 2 - 4s 2 X 2 , 2 (x, = (f + 2sx) 2 - 4s 2 X 2 . (4.22) 
Substitute (Olji . (1416]) into (T4T20D to get 

^3 = ( P " ^1*2^1^2) , (4.23) 

where 

P = 4s 2 (x 2 - x 2 )[2(r - J9 2 )x 2 - r 2 + r 4 ] + 8s 2 [(r - 2 X 2 )x 2 + t X 2 }£ 
- 2[(r - p 2 - 2s 2 )x 2 + r(p 2 - 2s 2 ) - r 4 ]£ 2 - 2r£ 3 - £ 4 . 

Denote 

Q = (f + r + 2s 2 - p 2 ) 2 - 4sV - (p 2 - 2s 2 ) 2 + r 4 . 
The following identity holds 

P 2 - <$> 1 ® 2 iZr 1 y 2 = 4x 2 (r£ 2 + ax 2 - rcx)Q 2 . 
Using the second equation (14.81) we find from (I4.23P 

w^—^iy^- yfP 2 ), (4.24) 

where 

P x = P + 2a;/xQ, P 2 = P - 2x/j5, /1 2 = /x^. (4.25) 
Simultaneously, from (14.81) . 

/i 2 = r £ 2 + ax 2 - ro-. (4.26) 

Thus, the equations (I4.15p . (I4.16p . (I4.18p . (I4.2ip and (I4.24p give the dependencies of all phase 
variables (12.151) on the two variables x, £ taken as the intermediate ones. 
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4.3. The change of variables. Introduce the local coordinates tx, t 2 according to (12.401) on 
the (x, £)-plane. The discriminant of the corresponding equation ( I2.39|) now coincides with 
the value /i 2 in ( I4.25P and is, therefore, non-negative. So the variables ti,t 2 are real. Solve 
(g3S} with respect to 

x h + t 2 h + t 2 ' ^ t x + t 2 

Here we denote for brevity 



U x = -a, U 2 = sjt\ - a. (4.28) 

Note that the signs of these radicals should be chosen arbitrary to supply, at given tx,t 2 , 
all possible triples x, in (I4.27P satisfying the system (I4.26p . (I2.40p . The following easily 
checked identities simplify some calculation below: 

t~\ — to , . 

xa=- -rf, 4.29 



(t x t 2 + a + UMihh + a- U X U 



2) 



(tl+t 



2) 



\2 



ht 2 + a- U X U 2 (Ux + u 2 



\2 



txh -a- U X U 2 (tx - t 



2 



\2 



a, (4.30) 
(4.31) 



We begin with fii, /12 defining the denominators in (I4.15p . (I4.16p . Put ip = 16s 2 r + \&i + 
ty 2 . According to (I4.14p we have 

A*i = t-(V>-2-s/*i*2) s A*2 = ^(^ + 2 v / ^i^ 2 ). 
8s 8s 

By definition ip 2 — 4\I' 1 \E' 2 = bAs 2 fi 2 , and we write 



Hence 



lbs 



lbs 

Substituting (I4.27P we get 
where 

(fx = V 2sy9 + E7i , <^2 = VZsy/r + U~ 2 , 

iPx = ^SyfF-Ux, 1p 2 = yJlSyfi-U*. 

This notation is intermediate, the signs chosen here do not affect the final expressions. The 
only important value is the new radical 



R= Vtit 2 + <T + UiU 2 . (4.32) 
Now the expressions for fix, ^2 become 

R 2 (ifm-^iip 2 ) 2 R 2 (ip 1 ip 2 + ^ 2 ) 2 

* = 4s(t 1 + t 2 ) 2 ' ^ = As(h+h) 2 • (433) 
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We can avoid "double roots" by using the expanded form: 

R 2 (As 2 r + U X U 2 - V X V 2 ) R 2 (As 2 r + U X U 2 + V X V 2 ) , . 

* = ' M2 = H^TW ■ (434) 

Here we introduce new algebraic radicals 



V 1 = J4s 2 x 2 -t 2 , V 2 = \Us 2 x 2 -t 2 (4.35) 



having arbitrary signs. 

To find x x ,x 2 , in addition to fl3~2ED, ( fP2]l . (1435]) denote 



M x = y/ti + t + r 2 , M 2 = Jt 2 + r + 



r 2 



X x = X + ^X 2 - 4r 2 $!$ 
X 2 = X - ^X 2 - 4x2$!$ 



(t x + £ 2 



2 • 

2 J1^2'Ar2 



ArR 2 M 2 N( 

2 



(tl+t 2 ) 2 ' 

For the numerators in (14. 15|) we obtain 

Ar\x 2 - r) + r( ± v 7 ^) 2 = ^(V^i ± V^f ■ 



Then from (Qijl . ( jUBJ , (14331 w e get 

/.U..V, • .U.A.A 

V ^1^2 + ^1^2 / 



2sr / .U, \ , • M X N 2 \ 2 2.s, / M 2 iVi - M X N 2 X 2 



= —5- , , , , x 2 



r 2 \ (p X (p 2 + i) x ip 2 J ~ r 2 \ ip x ip 2 - ip x ip 2 

or in the expanded form 

2st (t x + r){t 2 + r) - r 4 + M X N X M 2 N 2 



x x 
x 2 



r 2 4s 2 r + U X U 2 + V X V 2 

2st (t x + r)(t 2 + r) - r 4 - M X N X M 2 N 2 
r 2 As 2 t + U X U 2 - V X V 2 



(4.36) 



N x = V*i + t - ?" 2 , N 2 = ^jt 2 + T- r 2 . 

The signs here are also arbitrary. For the polynomials ( 14. 1 7|) we have 

2R 2 M 2 M 2 ^ 2R 2 N 2 N 2 

Let X = Ar 4 (x 2 - r) + r($i + $ 2 ). Then from (14371) and (14301) we find 

„ 2jR 16r 2 r 4 (t! - t 2 ) 2 i? 4 

X 2 - 4r 2 $ 1 $ 2 = — . 

(ti + 1 2 ) 4 

Hence 2r A /$i$ 2 = ^ X X X 2 , where 

AtR 2 N 2 M 2 



(4.38) 



(4.39) 



The expressions for y x ,y 2 in terms of t x , t 2 are now easily found from the definition (14. 6[) 
of n x ,n 2 . Write 



= r 2 xi;U 2 - /i 2 , ry 2 A*i = r 2 x 2 /ii - fx 2 . 
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Then after substitution of lOlj) . (jCTj) , ( fQ7|» we get 

rgi + t 2 - 2p 2 + 2r) - ^ + M^M^ 
Vl ~ S As 2 t + U x XJ 2 + V x V 2 ' f44 s 

_ rfti +t 2 - 2p 2 + 2t) - U X U 2 - y - } 

V2 ~ S 4s 2 r + U t U 2 - V 1 V 2 

The variables Zi,z 2 are expressed from (14.171) . (I4.18p : 

R \I\ A I-> + N X N 2 R M X M 2 - N X N 2 

Zl ~ ~JTr h + t 2 ' Z2 ~ ~JT r h + h ' (4 ' 41) 

To find the components of the angular velocity, start with «7 3 . Using the identity (14.301) 
represent the polynomials P, Q as follows 

4£ 2 ~ 2£ ~ 

P = i — ^^P, Q = — —Q- 

Here 

P = -(t 2 - r 4 )(t 2 - r 4 ) + r[(2s 2 - t 2 )t\ + (2s 2 - t x )t\ - p 2 (t 2 + t\) 
+r\h +t 2 - As 2 + 2p 2 )} + 2r 2 (2s 2 - p 2 )^ + t 2 ) 
+T \t l +t 2 + As 2 -2p 2 ) + r\ 

Q = t x t 2 + (2s 2 - p 2 )(t x + t 2 ) + r 4 + rfa + t 2 + As 2 - 2p 2 ) + r 2 . 

In virtue of (I4.35[) . (14.361) we obtain for the functions (I4.25P 

Pl = wf^ MlNlvl A = wf^ M ^ v * ■ 

Then from (14. 24ft using the identity (I4.29P we find 

lh - U 2 M 2 N 2 Vi - MiN^ 1 M 2 N 2 Vi - MiNiV 2 

W3 ~ V2^ t\-t 2 U 1 + U 2 ■ {AA2} 

Note that the expressions for Xi, yi, Z{ contain only the product V\V 2 . Therefore we can write 
the signs of the terms in the numerator of W3 in different ways. Then, having fixed the form 
(14.421) and applying (14.211) to find w x ,w 2 , we must point out the rule to define the signs of 
the radicals \/0i, \/® 2 in a way that guarantees all combinations satisfying the system of 
three equations ( I4.2p linear in Wi, together with ( I4.42p . Since the determinant of this system 
is zero by virtue of ( 14. 3p . we need to check only two of these equations, e.g., 

(y 2 + 2s)wi + x\w 2 + Z1W3 = 0, 
x 2 wi + (yi + 2s)u> 2 + z 2 w 3 = 0. 

From (I4.38P we write 

V2st M 2 N l + M x N 2 ^_ V2st M 2 N l - M X N 2 

X l = , , , > V X 2 



(4.43) 



r ip x ip 2 + ^1^2 " r ip!(p 2 - ^1^2 

Here the signs are chosen to satisfy the necessary condition ^Jxl^fx^ = x, where the value x 
is defined according to (14.271) . For 61, 02 the equations (I4.22p and (I4.27P then give 

(h + h) 2 ^ = -2R 2 <p 2 ipi (h + t 2 ) 2 e 2 = -2R 2 ^ 2 V 2 . 

Hence, introducing £i )2 = ±1 we write 

(ti + t 2 ) Je[ = i £1 y/2R^ 2 , (tx + t 2 ) JQ~ 2 = i e 2 s/2R^ 2 . 
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Then 

rR (e 2 (fl - e^g)Vi + (e^f - £ 2 ^)^> 
Wl ~ As^f (h + t 2 )(MiJV 2 - M 2 N 1 ) 

rR (e 2 <pl - e^gyVi - (giyj - £2^1)^2 

W2 ~ As^ (ti + t 2 )(MxN 2 + M 2 N 1 ) 

Choose £2 = — £i- The expressions 

rR V1-V2 rR V 1 + V 2 

wi = ±— FT , — —r^mrr tttt. w 2 = =F- 



■ v^(ti + £2) M 2 JVi - MiiVa ' v^(ti + t 2 ) M 2 Ni + M X N 2 ' 

do not satisfy (I4.43P for any choice of the signs. For E\ = e 2 = e out of two possibilities 
e = ±1 only one satisfies (I4.42p and (I4.43p . namely, 

= r(U 1 V 2 + U^R 

Wl ~ 2 Sv ^(t! + t a )(M 2 JV 1 - M X N 2 ) ' u m 

r{U x V 2 - U 2 V X )R ^ A ^> 

w 2 



28y/sr{t 1 + t 2 )(M 2 A^ 1 + M X N 2 ) ' 

To accomplish the algebraic solution note that we still have a "double root" expression 
for R. In fact, it can be eliminated. Put 



X = y/O. 

This value can be either real or pure imaginary. Denote 



K x = Jh + x, K 2 = Jt 2 + x, Li = Jti - x, L 2 = ^t 2 - x. 

Then 

R=^ 7 ={K 1 K 2 + L 1 L 2 ), U X = K X L X , U 2 = K 2 L 2 . 



Finally, the formulae (jCTjl . (OOjl . (jOTp . (jQ2|l and flCTj) algebraically define the 
values of complex phase coordinates (12.151) for given ti,t 2 up to the choice of signs for the 
following radicals 

V^, K u K 2 , L u L 2 , V u V 2 , M l5 M 2 , iV 1; JV 2 . 

Here V^, Mj, can be real or pure imaginary, and the pairs K\, K 2 and L\, L 2 can be real or 
complex conjugate. Accessible regions for t\,t 2 are defined by the condition that the initial 
phase variables aj,(3j,Uj (j = 1,2,3) in (12.151) are real. 

4.4. The equations of motion and the real solution. To derive the differential equa- 
tions, which describe the time evolution of the variables ti,t 2 let us use the variables s 2 
as the intermediate ones. From (I2.35p . (14 .7p we get 



51 = —2^' 52 = —2^- (445) 
Calculate the derivatives by virtue of (13. 3ft . Then from (13. 4ft we obtain 



ds\ 

~dt = i 4x^^ Zl + Z2 ^ XlW2 ~ X2Wl ^' 
~dt = *4^3 ^ ~ Z2 ^ XlW2 + X2Wl) - 



(4.46) 



On the other hand, 



dsi dsidti dsidt 2 ds 2 ds 2 dti ds 2 dt 2 . . 

~dt ~ ~dt[~dt ' + ~dt^~dt' ~dt ~ dt^H + ~dT 2 ~dE' ^' ' 
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Substitute fl427j) to f l4~45j) to obtain 

ds x t x t 2 - o + U X U 2 Ml ds x _ t x t 2 - a + U X U 2 Ml 

dh = ~ 2^(t x -t 2 ) 2 Th' dh~ 2^F(t 1 -t 2 ) 2 

ds 2 tjh -a + U X U 2 Nl ds 2 _ t x t 2 - a + U X U 2 Nl 

dt x = 2^{t x -t 2 ) 2 fV dt 2 ~ 2 v /7(t 1 -t 2 ) 2 U 2 

From (1PTD . (OS]) we have 

d*l = y/r(ti -t 2 )?7i / M 2 rfg 2 _ jyajglN 

dt r 2 (t x t 2 - a + U x U 2 y 1 dt 1 dt J ' 

C^2 y/r (tj - t 2 ) U 2 2 ds 2 2 dSj 

dt r 2 (tit 2 -a + C/i[/ 2 ) 1 2 dt 2 dt h 
Express the values flOB]) in terms of ti,t 2 using ffl~39]) . f H~4Tj) . fl4~44j) : 

dsi _ . (t x t 2 + tr + U 1 U 2 )M X M 2 (M 1 N 2 V 2 - M 2 N X V X ) 



(4.48) 



(4.49) 



dt 2^r(U x -U 2 ) 2 (t x -t 2 ) 
ds 2 . (t x t 2 + a + U 1 U 2 )N 1 N 2 (M 2 N 1 V 2 - M X N 2 V X ) 



dt 2y/2sT{U x -U 2 ) 2 {t x -t 2 ) 

Substitute these expressions to (14.491) . Then in view of the relation f!4.3ip after some obvious 
transformations we come to the equations of the Kowalevski type 

(*i " *2)^ = \j^s 2 x 2 -mi-a)[^-{t l + TY] , 

(h - = \/^( 4s v - mi - * 4 - (* 2 + rf] . 

To obtain the expressions for the real phase variables use the inverse linear mapping to 
the change fl235) . From (TOHD, (jODjl . (T4T4TD we find a jy (j = 1, 2, 3) as functions of t x , t 2 : 



- r 2 U x U 2 )(As 2 r + U X U 2 ) - (r + / - ) .U : A ; .U_> .V L ,\ , \ 2 
4r 2 sr(f/! + ?7 2 ) 2 

. (4 - r 2 ^iZ7 2 )ViV2 - (4s 2 r + U x U 2 )(t + r 2 )M x N x M 2 N 2 
a 2 = i 



tt'3 



Ar 2 sr(U x + U 2 ) 2 
R M X M 2 



, y 2h + t 2 

, _ . (B + r 2 LW 2 ) W 2 - (4s 2 r + C/iC/ 2 )(r - r 2 )M x N x M 2 N 2 { ^ 0U) 
Pi — 1 



2 



Ar 2 s t(U x + U 2 ) 2 
(B + r 2 U x U 2 )(As 2 T + U X U 2 ) - (r - r 2 )M 1 AT 1 M 2 AT 2 V r 1 \/ 2 



o . i2 N X N 2 

P3 = -1- 



4r 2 sr([A + U 2 ) 



r^/2 t x + t 2 

Here for brevity we put 

A = [(t x + r + r 2 )(t 2 + t + r 2 ) - 2(p 2 + r 2 )r 2 ]r, 
fi = \(t x + r-r 2 )(t 2 + T- r 2 ) + 2(p 2 - r 2 )r 2 ]r. 
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The angular velocities Uj (j = 1,2,3) are found from (12TT5D . flQ2|) . fl444jl : 

1 Ars a/st if — t 2 



iR M 2 N l U 2 V 1 + M l N 2 U 1 V 2 , . 

W2 = -W^ tFfl ' (4 ' 51) 

t/i - £/ 2 M 2 iV 2 Vi - MxiViVs 

These explicit expressions of the phase variables in terms of the t±, t 2 along with the separated 
differential equations give the complete analytical solution for the third critical subsystem. 



5. Conclusion 

In this article we presented the results dealing with analytical solutions in the generalized 
Kowalevski problem. The questions of qualitative and topological analysis of the arising 
subsystems are connected with the investigation of the multi-valued dependencies (13.14p 
and (14.501) . (I4.5ip . Formally these dependencies are the inverse mappings for the coverings 
of degree 2 n , where n is the number of radicals with arbitrary signs. Lately, new methods 
of investigating the phase topology of algebraically solvable systems were proposed in [30] 
based on calculating the invariants of some Z 2 -linear mappings. (When this article was 
already published a theory of the topological analysis based on Boolean vector-functions 
was presented in [JTJ H2] with applications to the classical problems and to the systems 
described above.) 

Up to this moment the first critical subsystem (the Bogoyavlensky case) has not received 
any analytical solution despite the fact that its topology is completely investigated in [15] 
and various Lax representations are found in [8[ [13] . For the case of the Kowalevski gyrostat 
in two constant fields no analytical solutions are found except for the special periodic motions 
and their bifurcations studied in [32] . though all critical subsystems are analytically described 
in [23]. 

The Kowalevski problem in the rigid body dynamics remains as challenging as it was 
during the last 120 years. 

The Russian version of this article was written with the financial support of the RFBR 
and the Volgograd Region Administration (grant No 10-01-97001). 
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